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1. Introduction and Results. 

In this letter, we use the methods recently developed in [21 El to prove 
that the integrated density of states (IDS) Nx{E) for a random Schrodinger 
operator H^{X) = Hq + XVuj is a uniformly Holder continuous function of the 
disorder parameter A at energies E for which the unperturbed operator Hq 
has a continuous IDS Nq{E), under fairly general conditions. Moreover, the 
uniformity in A implies that N\{E) — Nq{E) is (Holder) continuous in A, as 
A — > 0, at points E of (Holder) continuity of Nq{E). This result applies to 
random Schrodinger operators on the lattice Z'^ and on the continuum IR'^, 
given as perturbations of a deterministic, background operator Hq = (— iV — 
^o)^ + ^0- We assume that the background operator is self-adjoint with 
operator core C^{X) (smooth, compactly supported functions on X) for 
X = Z,'^ OT X = . For simplicity, we assume that > —Mq > — oo, for 
some finite constant Mq. In addition, we require that Hq is gauge invariant 
under translations by elements of Z"^. Specifically, this means that for every 
m G Z"^, we have Vo{x + m) = Vb(x) and Aq{x + m) = ^0(2;) + V(j)m{x) 
for some function (prn- For X = Z'^, the operator (— — ^0)^ represents a 
short-range, e.g. nearest neighbor, hopping matrix. 

We consider Anderson-type random potentials constructed from a fam- 
ily of independent, identically distributed (iid) random variables {uij \ j G 
Z'^}. On the lattice the potential acts as 

(VL/)(m) = iOmfim), m G / G (1.1) 

On IR'^, the potential Vuj{x) also depends on the single-site potential u, and 
is also a multiplication operator given by 

iV^f){x) = ^A^ - j)fi^)^ f e L\1R''). (1.2) 

Precise hypotheses on the single-site potential u and the random variables 
{ujj I j G Z'^} are given below. 

The family of random Schrodinger operators is given by 

HUX) = Ho + XV^. (1.3) 

The parameter A > is a measure of the disorder strength, and we consider 
the other parameters entering into the construction of V^, that is ||ti||oo 
and the distribution of loq, as fixed. As we are interested in the explicit 
dependence on A, we will write Hx for Hi^{X) and suppress lo in the notation. 
Due to the assumed gauge invariance under shifts of Hq, and the explicit 
form of the random potential given in (|1.1|) and (|1.2j) . the random operator 
H^{X), for fixed A, is ergodic with respect to the gauge twisted shifts 

V'(x) e^'^'"('^V(a; - m) , V € ^^(X). (1.4) 

We mention that the results of this note are easily modified to apply 
to the random operators describing acoustic and electromagnetic waves in 
disordered media, and we refer the reader to El El • 
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Our result follows the investigation initiated in 2^ where a proof of the 
Holder continuity in energy of the IDS is given that relies on the continuity 
of the IDS for the unperturbed, background operator Hq. As in the first 
part of |2j, we require that the IDS Nq{E) for the background operator Hq 
exists and that it is Holder continuous in the energy. The proof is local in the 
energy and applies at any energy E at which No(E) is Holder continuous. In 
particular, it applies to Landau Hamiltonians away from the Landau levels, 
where No(E) is discontinuous. 

Before stating our results, let us make precise the hypotheses on the 
random potential. 

Hypothesis (HI). The family of iid random variables {ujj \ j € Z"^} is 
distributed with a density h € L°°{IR) with compact support. 

Hypothesis (H2). The single-site potential u >0 is bounded with compact 
support. There exists an open subset O C supp u and a positive constant 
n > so that > K > 0. 

We first recall a result of f2], that Holder continuity in energy of the IDS 
for Hq implies continuity of the IDS for Hi^(X), with a constant and Holder 
exponent independent of A. 

Theorem 1.1. We assume that the Schrddinger operator Hq admits an 
IDS Nq{E) that is Holder continuous on the interval I C IR with Holder 
exponent < gi < 1, that is 

\NQ{E)-NQ{E')\<CQ{qi,I)\E-E'\i\ (1.5) 

for all E,E' G /, and some finite constant < Co(gi,/) < oo. We assume 
hypotheses (HI) and (H2) on the random potential V^j. Then, for any con- 
stant < q < qiq* /{qi -\- 2), where q* = I for and < q* < I for 
L'^{IR'^) (see \2.1fffi ). there exists a finite positive constant C{q,I), indepen- 
dent of X, so that for any A 7^ 0, and any E, E' £ I, we have 

\N^{E) - Nx{E')\ < C{q,I)\E - E'\^. (1.6) 



Note that the exponent q obtained by this method is roughly 1/3 whereas 
it is believed that it should hold with q = 1 (see section 4). A similar result 
was obtained recently by one of us jl2j . using a method quite different from 
that in 0. 

We now present the main result of this note. 

Theorem 1.2. Under the hypotheses of Theorem 1.1, for any bounded 
interval J C M, there exists a finite, positive constant C{q,I,J), such that 
if A, A' € J, we have 



\NxiE)-Nx'{E)\<Ciq,I,J)\X-X' 



(1.7) 
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for any E £ I and < 52 < 2g/(g + 3), where < q < qiq*/{qi + 2) is the 
exponent in 

Until recently, it was not known that the IDS remained bounded in 
the weak disorder limit A — > 0. In particular, result (|1.7|) was known 
only for closed intervals J disjoint from zero. This result follows from the 
Helffer-Sjostrand formula (see section 3 and also However, the constant 
C{qi,I,J) obtained from that proof scales like [dist (J, 0)]~^. 

Recall that control of the IDS comes from the Wegner estimate, 

iP{dist {a{HA),E) <7]}< Cg{X)\A\rj'^, (1.8) 

for any < g < 1. Here H\ is the restriction, with suitable boundary 
conditions, of -ff(^(A) to a bounded open set A of volume |A|. In the usual 
proof of the Wegner estimate the constant Cq{X) diverges as 1/A as 

A ^ 0. In [2], a different proof of the Wegner estimate is given for which the 
constant C{q,I) is uniformly bounded in A. The only deficit of this proof 
is that the Holder exponent q for the IDS Nx{E) must be taken sufficiently 
small (as stated in Theorem 1.1) relative to the assumed Holder exponent 
< (?i < 1 of the IDS of Nq(E) in (|1.5)) . In particular, the bound gives 
no information about the density of states (DOS) p\{E) = dN\{E) / dE (see 
section 4 for a further discussion of the DOS). 

We have the following two corollaries of Theorems 1.1 and 1.2. 

Corollary 1.3. Under the same assumptions as Theorem 1.1, let J C M 
be any closed, bounded interval containing 0. Then, there exists a finite, 
positive constant C{q,I, J), so that we have for any E E I and A G J, 

\Nx{E)-No{E)\<C{qi,I,J)\X\'^\ (1.9) 

where < 52 < 2g/(g + 3), where < q < qiq* /{qi + 2) is the exponent in 

There is a version of Theorem 1.2 and Corollary 1.3 with the weaker hy- 
pothesis of continuity for Nq{E) and with a correspondingly weaker result. 

Corollary 1.4. We assume that the Schrodinger operator Hq admits an 
IDS Nq{E) that is continuous at E. Then, under the same hypotheses (HI) 
and (H2) as in Theorem 1.1, we have for any A that the IDS N\ is also 
continuous at E and that 

lim Nx'{E) = Nx{E). (1.10) 

In general, as the IDS Nq{E) is a monotone increasing function, this result 
applies at all but a countable set of energies. 

In section 2, we recall the proof of Theorem 1.1. The proofs of Theorem 
1.2, and Corollaries 1.3-1.4, are given in section 3. We conclude with some 
comments about the behavior of the density of states in section 4. While 
preparing this letter, we learned that Germinet and Klein ^ have proved 
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a version of (|1.7)) for intervals J away from zero. We thank F. Germinet 
(private communication) for showing us the use of (|.S.9j) that improves our 
original estimates on q2- 

It is clear that there are various generalizations of our results. For exam- 
ple, hypothesis (HI) can be weakened to allow unbounded random variables 
with the first two moments bounded. 

2. Sketch of the Proof of Theorem 1.1. 

For completeness, let us sketch the proof of Theorem 1.1 that appears in 
[21. We assume hypotheses (H1)-(H2) and condition ((T3|) on the IDS Nq{E) 
for the background operator Hq. Let A C / be a sufficiently small closed 
interval, and let A D A be a bounded interval with |A| = OdAj"), for some 
a G (0,1). First, one proves that ()1.5|) implies that for all A sufficiently 
large, depending on A, there exists a finite constant Ci(/,d) > so that 

TrE^{K) < Ci{I,d)\A\'''\A\. (2.1) 

Next, we consider the local spectral projector E\{A) for H\ and write 

TrE^iA) = TrEA{A)E(^{A) + TrEA{A)E^{A^), (2.2) 

where A'^ = IR\A. The first term on the right in (|2.2j) is easily seen to be 
bounded by 

Tr Ea(A)S^(A) < TrE^{A) < C7i(/, d)|A|°«i |A|, (2.3) 

and is already of order |Aj^, for any q < aqi. 

The second term on the right of (|2.2() is estimated in second-order per- 
turbation theory. Let £^ E A be the center of the interval A, and write 

TrEA{A)E^{A'^) = TrEA{A){HA - E)E^{A'^){H^ - E)-' 

-XTrEA{A)VAE^{A'){H^ - E)-' 

= i^ + iii)- (2.4) 

Since the distance from A'^ to E is of order |A[", we easily see that term (i) 
of (|2.4j) is bounded as 

|(i)| < |A|i-"rri?A(A), (2.5) 

so that as < a < 1 and |A| < 1, we can move this term to the left in 1)2.2(1 . 
Continuing with (ii), we repeat the calculation in ()2.4j) . now to the left of 
E\{A), and obtain 

{ii) = -XTr{HA- E)Ea{A)VaE^{A^){H^ - E)-^ 
+\^TrVAEA{A)VAE^{A-){H^ - Ey^ 
= (m) + (iv). (2.6) 

Term (in) is estimated as in (|2.5() and we obtain 

|(m)| < A|A|i-2"||T4||Tri?A(A), (2.7) 
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where Vx is the potential obtained by replacing ujj by the maximal value of 
\ujj\. Term (iv) in ()2.6|) can be bounded above by 

\{iv)\ < X^\A\-^''TrVAE^{A)VA. (2.8) 

Taking the expectation and replacing by the upper bound V^, we find 
that we must estimate 

IE{Tr{VlEA{A))}. (2.9) 

This is done using estimates on the spectral shift function comparing the 
two local Hamiltonians with one random variable fixed respectively at its 
maximum and minimum values. For the lattice case, this is a rank one 
perturbation, so the corresponding spectral shift is bounded by one, the 
rank of the perturbation. For the continuous case, the perturbation is no 
longer finite rank, but we may use the local L^-estimate on the spectral shift 
function proved in 3 . In either case we obtain 

IE{Tr{VlEA{A))} < C^{I,q* ,u)X-'\Af \A\, (2.10) 

where the exponent q* in (|2.1fl|) is i) q* = 1 in the lattice case, ii) < < 1 
in the continuum. 

As a consequence, term (iv) in (|2.6|) can be bounded by 

lE{\{iv)\} < X\Af-^''Ci{I,d)C4{I,q*,u)\A\. (2.11) 

Putting together (|^ . 1^7]) . and (EHU, we obtain 

{1 - |A|i-" - X\A\'^-^''\\VA\\}lE{TrEA{A)} 

< {X\Af~^''Ci{u,d)C4{q*,u) + C2\Ar')\A\. (2.12) 

By choosing the optimal a < 1/2, it is clear from this expression that 
Theorem 1.1 holds with 0<q< qiq*/{qi + 2). □ 

3. Proof of Theorem 1.2. 

The almost-sure existence of the IDS for random Hamiltonians of the type 
considered here is well-known and we refer the reader to ^ \U)\ lllj . The 
IDS N\{E) is given in terms of the spectral projector P\{E) associated with 
Hi_j{X) and the interval {—oo,E] C IR. For the lattice case, with Hilbert 
space the IDS Nx{E) is given by 

Nx{E) = lE{Tr6oPxiE)6o} 

= IEmPx{E)\0)}, (3.1) 

where is the function supported at and \x) is the state at site x € Z"^. 
For the continuous case on IR'^, the IDS Nq{E) is given by 

Nx{E) = IE{TrxoPx{E)xo}, (3.2) 

with xo the characteristic function on the unit cube in IR'^. To unify the 
notation, we will write xo for the characteristic function on the unit cube 
as in iS. in the continuous case, or for the projector 5q as in /TO) in the 
lattice case. 
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Proof of Theorem 1.2. Fix A, A' G J and ^ G /. Choose g G C^{IR), 
depending on E, A, and A', with < g < 1 and 

where < a < 1 wih be determined. The choice of g G obeying (|3.3j) is 
basicahy arbitrary, however we require that 

||ff(^')||oo<C7|A-A'|-^", (j = 1,2,4), (3.4) 

with some constant independent of E, A, and A' (this can be done). 
We have 

Nx{E)-Ny{E) = Nx{E) - IE {Trxo[g{Hx)]\o} 

+IE {Trxo[9{Hx')?Xo} - Ny{E) 

+IE {TrxoMHx)? - [g{Hx')f)xo} ■ (3.5) 

The monotonicity of Nx{E) with respect to energy and the properties of g 
imply that 

IE {Trxo[g{Hx)?Xo} Xo < NxiE + |A - A']")- (3.6) 

It follows from Theorem 1.1 that 

IE {Trxo[g{Hx)]\o} - Nx{E) < Nx{E + \X - X'n - Nx{E) 

< Ciqi,I,J)\X-X'r, (3.7) 

for any < q < qiq*/{qi + 2), and an identical estimate holds for the second 
term on the right in 1)3. 5() . 

It remains to estimate the last term on the right in 1)3. 5|) . Using the 
identity 2{A^ - B^) = {A{A -B) + {A- B)A} + {B{A -B) + {A- B)B}, 
we can write the last term in ()3.5|) as 

IE{TrxoMHx)?-[giHy)?)Xo} 

= IE {Trxog{Hx){g{Hx) - g{Hy))xo} 

+ IE {Trxog(.Hy ){g{Hx) - g{Hx'))xo} , (3-8) 

where, to reduce the number of terms, we have made use of the following 
identity: If A^^ and B^ are .^"^-ergodic operators such that XoA^jB^^xo is 
trace class, then we have 

IE {TrxoA^B^xo} = IE {TrxoB^A^xo} ■ (3.9) 

(We use this identity in a more crucial way below.) We note that the trace 
norm ||xo9(-f^A) |ji is bounded uniformly in A G J as well as in the random 
couplings ojj. In the continuum, we are using here that Hx is bounded from 
below. 

We express the difference {g{Hx)—g{Hx')) in terms of the resolvents using 
the Helffer-Sjostrand formula, which we now recall (see, for example, ^.6, for 
details). Given / G CQ^^(iR), we denote by fk an almost analytic extension 



of / of order k, which is a function defined in a complex neighborhood of 
the support of / having the property that + iO) = f{x) and that 

\d,fk{x + iy)\ ~ \f^'+'\x)\\yt as |y| ^ 0, (3.10) 

where dz = dx + idy. For the construction of such a function, which is not 
unique, we refer to j^. Let R\{z) = [Hx — z)~^ denote the resolvent of H\. 
For functions g as in (|3.3() . the functional calculus gives 

5(Fa) - 5(^A') = / d,g{z)Rx{z)V^Rx'{z)(fz, (3.11) 

71" Ja 

with ^ an extension of order 3 (recall that g E C^). 

Let us estimate the first term on the right in ()3.8() . The estimate for the 
second term is similar. We substitute the Helffer-Sjostrand formula 1)3. 11() 
and find 

Trxog{Hx){g{Hx) - g{Hx'))xo 
_ (A - A') 
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dz~g{z)Trxog{Hx)Rx{z)V^Rx'{z)xod'z. (3.12) 



Using the second resolvent identity, we rewrite the operator involving resol- 
vents as 

Rx{z)Vu,Rx>{z) = Rx{z)Vu,{Rx'{z) - Rx{z)) + Rx{z)Vu,Rx{z) 

= {\' -X)Rx{z)V^Rx'{z)VMz) 

+Rx{z)V^Rx{z) 
= {i) + {n). (3.13) 

The integral in (|3.12j) involving the first term [i) in (|3.13j) is estimated 
as follows. The resolvents are bounded by |Imz|~^ as |Imz| — > 0, but this 
divergence is canceled by the estimate (|3.10|) for dzg (since we take an ex- 
tension of order 4) . Recalling the estimate (|3.4|) on the derivatives of g and 
noting |supp g'\ ^ 6^ , with = |A — A'|", we find that 

\d2g{z)\\lmz\-^(fz < C\X - A'|^3", (3.14) 



'a 



since we obtain a factor of 5^ ^ from the derivatives and a factor of 5i from 
the size of the domain of integration. Consequently, we find 



IE 



^^^^ / dz~g{z)Trxog{Hx)Rx{z)V^Rx'{z)VMz) 
71" Ja 

< CIE{\\xo9{Hx)\\i\K\L} |A - Af < c,\X - Af "3". (3.15) 

To evaluate the integral involving the second term (ii) of (|3.13|) . we apply 
H3.9() to the operator integrand in (ii) of 1)3. 13() . Inserting this into (|3.12|1 . 
we obtain for the integrand 

IE {Trxo9{Hx){Rx{z)VMz)xo)} = IE {Trxo9{H,)R^{zfV^xo} ■ (3.16) 



The integral becomes 

[ d^g{z)Trxog{Hx)Rx{zfV^Xod^z = -Trxog' {Hx)V^Xo- (3.17) 
Ja 

As a result, we obtain the following estimate for the term involving (ii) 

\lE{Trxog'{Hx)V^Xo}\ < Ci\X - AT"- (3.18) 

Combining the estimates (|3.7() . ()3.15() . and ()3.18|) . we obtain the upper 
bound for the right side of (|3.5() . 

\Nx{E)-Ny{E)\ 

< 2C{qi,I)\X - AT"^ + Co\X - A'|2-3" + Ci\X- Af (3.19) 

Comparing the exponents of |A — A'| in H3.19() . we can take < a < 1 so 
that aq = 2 — 3a, giving the exponent 2q/{q + 3). □ 

The proof of Corollary 1.3 follows simply by taking A' = 0. The continuity 
result of Corollary 1.4 is proved as follows. 

Proof of Corollary 1.4. It suffices to note that the proof of Theorem 1.1 
in |2] can be extended to prove that if Nq{E) continuous at E, then so is 
Nx{E). To see this, fix € iR at which Nq is continuous. Following the 
argument of [2], we see that the finite- volume estimate ()2.1() becomes the 
following. For any e > 0, there exists r] > such that for A = [E — r], E + rj], 
and all A sufficiently large, one has 

TrE^{A) < e|A|. (3.20) 

Without loss, we assume that < e, since the left side of (13.20(1 is non- 
increasing in 7]. Choose a closed sub interval A = [E — 7]^ , E + r]P], with 
p > 1. Following the argument in section 2 with this choice of A and A, the 
estimates (|2.5() . (|2.7() . and (|2.11l) now have the form: 

I(i)| <7?*^-irri?A(A), (3.21) 
< CoXr]P-^TrEA{A) (3.22) 
lE{\{iv)\} < CiXriP'i'-^\A\, (3.23) 

where the constants Co and Ci are independent of A, r/, and e and the 
exponent q* appears in (|2.1()|) . These imply that (j2.12|) has the form 

{1 - rf~^ - CoXr]P'^}IE{TrEA{A)} < CiXif^'-^ + e)|A|. (3.24) 

If we pick p > 3/(7*, then for sufficiently small e, we get for all A sufficiently 
large 

IE{TrEK{A)} < eC\A\, (3.25) 

for some finite constant C > since rj < e. This shows that the IDS Nx is 
continuous at E. 
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To complete the proof of Corollary 1.4, we return to equations (|(-i.5j) and 
(|3.7j) . We use the continuity of Nx(E) to control the first and the last terms 
on the right in (|3.5j) . For example, we need to estimate 

Nx{E) - IE {Trxog{Hx)xo} • (3-26) 

The monotonicity of N\{E) with respect to energy, and the properties of g, 
imply that 

IE {Trxog{Hx)xo} < Nx{E + |A - AT )• (3.27) 
It follows from the continuity and monotonicity in E of Nx{E) that 

< IE {Trxog{Hx)xo} - Nx{E) < Nx{E + |A - A'D - Nx{E) 0. 

(3.28) 

The estimate for the middle term of (|3.5|1 remains the same. Consequently, 
we have that 

lim Ny{E) = NxiE), (3.29) 

A' — >A 

at any point E of continuity of A'o, proving Corollary 1.4. □ 

Remark: This proof shows that in general one can control the modulus of 
continuity for the IDS Nx{E) of the random model using that of the free 
model. 



4. Additional Comments and Conjectures 

In certain situations, we are able to obtain more information about the 
density of states (DOS) px{E) and its behavior as A ^ 0. The DOS is the 
derivative of the IDS Nx{E) with respect to energy. Since the spectral shift 
function is pointwise bounded for the lattice model, it follows from j2j that 
the DOS is bounded except at possibly a countable set of energies. In this 
case, the DOS is given by 

Px{E) ^ ^(E) = l\mIE{lm{0\{Hx -E- .e)-i|0)}. (4.1) 

Let us suppose that the random variables ujj are Gaussian with mean zero. 
In this case, the almost-sure spectrum of -ff(^(A) is IR, for A 7^ 0, and the 
spectrum of Hq = A is [—2d, 2d]. If E G IR\[—2d,2d], the resolvent can be 
expanded in a Neumann series, 

00 

RxiE + ie) = Y^ Ro{E + ie) [-XVME + ie)f . (4.2) 

k=0 

The matrix elements of the free resolvent decay exponentially by the Combes- 
Thomas argument. Let do{E) be the distance from the spectrum of Hq to 
E. We then have the bound, 

\{x\Ro{E)\y)\<^^e-'om^'y\/2_ (4.3) 
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We take the expectation of the zero-zero matrix element in (|4.2j) . We expand 
the potentials and use the estimate (|4.3j) to control the sum over sites. 
We easily see that the power series converges absolutely provided 

\X\lE{\u;o\}Ci{d) 

where the constant Ci{d) depends on Co in 1)4. 3|) and the dimension. For 
example, for all A < 1, we have the convergent expansion 

oo 

px{E) = X'p^^He) + Y,^'p^''He), (4.5) 

fc=3 

for all 1^1 > [Ci{d)IE{\uJo\}]'^/(''+'^^ + 2d. 

This result, and the results on the IDS in this note, are steps towards 
proving the general conjecture concerning the regularity of the DOS. In par- 
ticular, under the hypotheses (II1)-(II2), we expect that the IDS is Lipschitz 
continuous, that is, we have g = 1 in H1.6|) . with a constant independent of 
A. Furthermore, if the unperturbed operator Hq has a Lipschitz continuous 
IDS, then we expect that 

\px{E)-Po{E)\<C,\X\^, (4.6) 

for some constant < Cg < oo, independent of A, and some < q < 1. 
Finally, if the distribution function for the random variable uiq is sufficiently 
regular, we expect that the IDS is also regular. 
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